Spaceborne NO2 observations are sensitive to coal mining and processing in the largest coal basin of Russia

Coal use exacerbates several major environmental problems including build-up of greenhouse gases and air quality deterioration. Although Kuzbass (Siberia) is one of the largest exploited coal basins worldwide, the role of regional coal mining and processing in atmospheric pollution is unknown. We outlined the Kuzbass coal basin by spaceborne night-lights and revealed a regional, long-term tropospheric NO2 anomaly (2005–2018) by spaceborne NO2 column observations (hereafter ‒ NO2). The spatial agreement between NO2 and night-lights indicates that the anomaly is attributable to an agglomeration of coal quarries and the cities in Kuzbass, that are heavily reliant on coal. A positive relationship between NO2 and interannual coal production suggested that the anomaly was related to coal in Kuzbass; ~ 1.0% of annual coal production increase induced ~ 0.5–0.6% of NO2 enhancement. As coal production accelerated since 2010, NO2 exhibited strikingly similar annual increases over Kuzbass in 2010–2014 (7%) and 2015–2019 (15%), compared to 2005–2009. Conversely, Siberian cities lacking a coal industry followed the global trend of reducing NO2 for the same periods (−5% and −14%, respectively), driven by fuel combustion improvements. Overall, we demonstrated that coal mining, processing and utilization can induce distinct tropospheric NO2 anomalies, detectable from space.

The global reliance on coal has triggered infamously adverse environmental effects as coal has become the main contributor of the atmospheric greenhouse gas accumulation and the most atmosphere-polluting source of energy as well 1 . These adverse effects stem from direct emissions of coal mining (scraping and fracturing of coal from ground) as well as from indirect emissions from coal mining (fuel combustion of coal mining machinery), coal processing (production, conversion to coke, use in metal production) and coal transportation (emissions from coal-carrying vehicles) 2 . There is a growing corpus of studies, relying on spaceborne observations to monitor and constrain atmospheric emissions, originating from coal power plants [3][4][5] . However, only few investigations have considered spaceborne observations of atmospheric byproducts from areas of coal mining and processing. Previous studies have focused on developed countries such as the U.S. 6 and Australia 7 , but 84% of the world's coal mining and processing are in emerging and developing economies 8 , following the development classification of the International Monetary Fund 9 . Knowledge of the atmospheric byproducts of coal mining and processing is incomplete for these countries in comparison with developed countries. This is peculiarly undesirable, because officially reported emissions in these countries can be strikingly inaccurate [10][11][12] .
A growing research interest in Chinese coal point-sources 13 has provided previously unknown information about coal byproducts in the atmosphere. However, various other emerging or developing countries remain unexplored in this respect. The post-Soviet states (accounting for 7% of global coal production) are the most salient examples of such blind-spots, considering the presence of large coal-producing regions such as Donbass in Ukraine and, especially, Kuzbass in Russia. Kuzbass, located in southwest Siberia, is the largest coal basin in Russia and one of the largest in the world (~ 300 billion tons of accessible reserves), containing 33% of the world's known coal deposits 14 . Despite the global importance of Kuzbass, only Oparin et al. 15 have addressed the air pollution over the region by using spaceborne remote sensing of snow cover. Besides their indirect evidence of air pollution, detectable from snow cover in Kuzbass, no other empirical work has attempted to estimate the atmospheric pollution or composition over Kuzbass. The literature contains only fragmented information about the effects of coal mining on the atmosphere over this coal-rich area such as the indication about elevated  16 The authors have never addressed this local increase, but such hint about a potential NO 2 atmospheric enhancement over a major coal basin is intriguing because NO 2 is not a common indicator of a direct coal mining outgassing. Rather, it is attributed to indirect emissions of coal mining, stemming from fuel combustion, mainly by motor vehicles 17 but also by the heavy machinery and transportation vehicles involved in the mining, processing, and transportation of coal 18 . In this context, a regional study reported that in Kuzbass, huge amounts of coal are mined (by excavators) and transported (by haul trucks) using heavy machinery that relies on the inefficient combustion of diesel fuel 19 . Given the scales of coal mining in Kuzbass, such machinery generates massive amounts of NO x emissions that might increase tropospheric NO 2 to levels that are seemingly harmful to respiratory and cardiovascular systems 20 .
For this reason, our study elucidates a potential link between coal mining/processing activities, and atmospheric NO 2 over Kuzbass by using a set of spaceborne remote sensing observations: NO 2 tropospheric columns from the Ozone Monitoring Instrument (OMI), night lights from the Operational Line-scan System (OLS), and urban pixels from the Moderate-resolution Imaging Spectroradiometer (MODIS). To this end, we (a) investigated the statistical agreement between NO 2 and the cluster of anthropogenic activities in Kuzbass (the concentrations of urban areas, industrial clusters, and opencast coal mines in the region); (b) examined the sensitivity of NO 2 to coal production (from reported data) in Kuzbass; and (c) separately compared the long-term trends of NO 2 for Kuzbass cities and other cities in Siberia that have no links to the coal industry.

Results
Kuzbass is located in southwestern Siberia and lies within the Kemerovo administrative region of Russia (Fig. 1), which is specialized in mining and processing industries, occupying ~ 50% of regional economy. Kuzbass covers the area of ~ 26,000 km 14 , contains ~ 300 billion tons of coal 14 , and is responsible for 58% of coal produced in Russia and 70% of exported coal 21 . Most of the coal (63.8%) is mined in open pits 22 , and there are 90 such mines 23 . Compared with other coal-rich regions in Russia, Kuzbass is a hotspot of coal mining, processing and utilization because of its mild climate and near-surface coal seams, which are conductive for open-pit mining and are cost-effective 21 . The adverse environmental impacts of coal-related industry in Kuzbass are clear, including 2.5 billion tons of waste (50% of the solid waste in Russia), with 98% of this waste originated from mining activities 22 .
While the boundaries of Kemerovo region are administratively defined (Fig. 1), the boundaries of the Kuzbass basin can vary depending on the application. There is no robust, validated map of Kuzbass coal quarries, compiled by established surveying techniques. Given the small size of the study area and the visual prominence of quarries against the green vegetation of the surrounding taiga ecosystem, we used Google Earth imagery and manually marked the centers of the largest mines in the region (Fig. S2 23 ). These quarries were used to create a polygon by connecting all the marginal coal quarries within the geographical cluster in the center of Kemerovo region with most identified quarries (68/81) within (red line, Fig. 2). This cluster reflects the heartland of the regional coal basin and is hereafter referred to as Kuzbass.
As shown in Fig. 2, we validated the mapping of mines in Kuzbass by using the latest estimates of cloud-free night lights from OLS measurements, which are a proxy for human activity 24 . Areas with opencast mines operated during the night normally exhibit distinct local increases in the digital number (DN) of night lights 25 . We www.nature.com/scientificreports/ also used the map of urban areas from MODIS, based on the unique phenology of build-up areas 26 to distinguish cases of increased night lights originating from urban areas (when the lights align with the MODIS urban map) from cases of increased night lights from a coal mine outside a built-up zone (when there is no overlap of lights with the MODIS urban map). A lack of overlap of night-lights with the urban population data has been previously used to identify night light signals that are unrelated to cities 27 . Figure 2 demonstrates that the mapping of the Kuzbass quarries realistically reflected the regional patterns of urbanization and the coal quarry allocation. The combination of coal mine mapping with night lights and the MODIS urban map (Fig. 2) hints at three broad types of area modified by human activity in Kuzbass: (1) urban areas without mines, identified by strong night lights overlapping with the MODIS urban map. An example is Novokuznetsk, a major regional conurbation with an industrial economy, specializing in metal production 28 , but relatively weak underground coal mining. The second (2) is urban areas that include coal mines. The urban areas are identified as above, but there are also coal mines, identified within the outlined areas. Examples are Mezhdurechensk and Prokopyevsk, cities planned and evolved exclusively as coal mining centers. The third area type is (3) coal mining zones, characterized by moderate night light intensity and the presence of mines but with no overlap with MODIS urban pixels. Note As mentioned, a distinct increase of tropospheric NO 2 over southwestern Siberia had been evidenced in a previous study 16 , but remained uncommented by its authors. We calculated the average NO 2 for 2005-2018 and revealed two distinct NO 2 spatial enhancements (e.g., positive anomalies) over southwestern Siberia as well, namely, over Kemerovo region (Fig. 3a). These two major long-term NO 2 anomalies were centered approximately over the regional industrial centers of Novokuznetsk (the major southern anomaly; peak NO 2 ~ 4.25 × 10 15 molecule/cm 2 ) and Kemerovo City (the minor northern anomaly; peak NO 2 ~ 2.90 × 10 15 molecule/cm 2 ). The maximum values of both anomalies corresponded to large urban areas with a population of ~ 0.5 million people. Moreover, the spatial gradient of NO 2 anomaly, we identified ( Fig. 3a) is consistent with that reported by Jamali et al. 16 over this area. Specifically, the major southern NO 2 anomaly showed the greatest NO 2 concentration in the center, which gradually weakened toward the borders of the identified cluster of coal mines in Kuzbass (black outline in Fig. 3a). Tropospheric NO 2 over Kuzbass was somewhat high, with an average of 3.22 ± 0.52 × 10 15 molecule/cm 2 within the outlined cluster of coal mines.
Most interestingly, the surge in night lights ( Fig. 2) looks strikingly similar to the elevated NO 2 over the coal basin in Kuzbass (Fig. 3a). The borders of the elevated NO 2 did not align with the Novokuznetsk metropolitan area (white outline labeled Novokuznetsk in Fig. 3), but spatially coincided with the large cluster of coal mines. As the NO 2 southern anomaly was unevenly distributed over the cluster of coal mines, we estimated the quantitative agreement between this anomaly and the night-light intensity, which reflects local human activity. Figure 3b reveals a moderate agreement between night lights and NO 2 within Kuzbass coal basin (Pearson correlation coefficient, r = 0.63 at p-value < 0.005), thereby corroborating the association between human activity and the NO 2 anomaly in the region. Although the correlation between the intensity of night lights over Kuzbass and the NO 2 anomaly is imperfect (r = 0.63; Fig. 3b), this finding demonstrates that this cluster of mines is likely the main anthropogenic driver of NO 2 emissions in the region. Notably, previous studies that considered NO 2 from OMI and night lights 25,29 , did not find such a strong correlation.
To assess the spatial association between coal mining/processing, and the NO 2 anomaly over Kuzbass, we estimated the statistical agreement between them for 2006-2018. This period differs from that in the previous section due to the availability of coal production data. At first glance, there was no strong correlation between the annual coal production rates and annual averages of NO 2 over Kuzbass (r < 0.50). Although NO 2 did not exhibit any significant temporal trend in Kuzbass in 2006-2018 (p-value > 0.05 based on a Mann-Kendal trend test), coal production exhibited clear increasing trend in the same period (p-value < 0.05, slope = 5.99). As coal is not related to night lights (i.e., anthropogenic activity), but to the local topography. In fact, the eastern margin of the region has the highest elevation in the study area (600-1050 m above sea level, asl; Fig. S3 in the Supplementary Material). NO 2 generally accumulates under two topographic conditions: the hindered dispersion of pollutants due to low penetration of winds in a trough, and pollution being trapped near the ground surface by temperature inversions formed by surrounding mountains 44 . The center of Kuzbass (the strongest NO 2 anomaly, red in a) is subjected to these topographic conditions, as it is located in a trough (~ 250 m asl) surrounded by hills (Fig. S3). It is therefore likely susceptible to temperature inversions, but there were no regional meteorological studies to evaluate this hypothesis in detail. www.nature.com/scientificreports/ production is often a monotonically increasing quantity and it exhibited clear increasing trend in Kuzbass, the actual variability in production can be hidden in an inter-annual signal. Hence, we tested the inter-annual variability (IAV) of coal production, a method previously applied to similar monotonically increasing quantities such as fossil fuel emissions 30 . Fundamentally, IAV coal reflects the difference (%) between an annual estimate of coal production during year n minus the coal production of year n-1. We identified reasonable agreement (r = 0.60) between IAV coal and the annual average estimates of NO 2 over Kuzbass (Fig. 4) despite seemingly incomplete regional coverage by OMI observations during these years.

Scientific
Although the statistical correlation is apparent, it was weakened by a few prominent differences in trend between NO 2 and IAV coal , such as in 2007-2008 (Fig. 4), when the world economic crisis struck most countries including Russia. IAV coal fell below 0%, which is the sole decrease in coal production in Kuzbass over the 11 year study period. The decline in coal production was more rapid than the general weakening of socio-economic activity within urban areas. This phase difference was driven by the dramatic increase in the cost of transporting coal. In particular, railway transportation (99% of coal is transported from Kuzbass to other regions by train) reached 40% of the Russian coal price in 2008 31 . Moreover, declines in related industries (iron, steel, chemical, and power) abruptly constrained the supply of resources required for coal mining 31 . In contrast, activities in urban areas were not so immediately affected by these factors.
The importance of coal production as one of the major drivers of NO 2 anomaly over Kuzbass was further evaluated by analyzing statistical agreement between city-scale estimates of NO 2 and population count from the national inventories 32 . This evaluation was based on the knowledge that in most cities, population, not coal is the main driver of NO 2 33 . To this end, two types of cities were used in the analysis: the major cities of Kuzbass (black circles) and the neighboring major cities in Siberia (red circles) in Fig. 5. This analysis revealed two distinct patterns. First, NO 2 over most Kuzbass cities (> 2.5 × 10 15 molecule/cm 2 except Mezhdurechensk) was distinctly higher, compared with other Siberian cities (< 2.5 × 10 15 molecule/cm 2 except Novosibirsk; the largest city of Siberia). Second, the strong linear association between NO 2 and population was discerned for other Siberian cities (r = 0.83), thereby, confirming the common, population-related driver of NO 2 in these cities. Notably, the correlation between OMI-based NO 2 and population in these Siberian cities is higher than for cities in the U.S., Europe, China, and India 34 . However, despite the similar number of compared points, there was no such agreement for the Kuzbass cities (r = 0.31), indicating that population size was not the main driver of NO 2 in Kuzbass. Figure 5, illustrating the relationship between NO 2 and population in Siberia, indirectly suggests that the recorded NO 2 levels in Kuzbass would correspond to a city with a population of > 2.5 million inhabitants. Overall, the combination of such high NO 2 tropospheric concentration and the lack of relationship with population points to the existence of another driver of the NO 2 increase in Kuzbass, unrelated to population, which might be coalassociated activities including mining, processing, transportation, and utilization of coal.
Coal mining and processing activities were likely among the main drivers of the NO 2 anomaly over Kuzbass, given the direct evidence (the overlap between the NO 2 anomaly and the cluster of coal mines, and the correlation between coal production and tropospheric NO 2 over Kuzbass) and indirect evidence (tropospheric NO 2 over Kuzbass not being related to population). This is intriguing and raises the question of whether monotonically increasing coal production affected the strength of the NO 2 anomaly over Kuzbass. To answer this question, we relied on previous findings, which reported dramatic decreases of spaceborne-based NO 2 (40% for some cities in the U.S.) over urban areas worldwide in 2015, compared with 2005 due to technological improvements and stricter regulations of emissions 34 . As these regulations were mostly underpinned to vehicular and stationary emissions, we tested whether the Russian NO 2 hotspots were within this decreasing global trend of the NO 2 , by assuming that the coal mining/processing activities were not affected by these measures. If this hypothesis is correct, the NO 2 change over Kuzbass would differ from that over the Siberian urban areas outside Kuzbass (red  Table S5 Figure 6 shows the differences between each of the latter periods and the baseline period. Interestingly, two distinctly different patterns were identified. First, the Siberian urban NO 2 hotspots outside Kuzbass did follow the decreasing global trend as NO 2 in 2010-2014 has been reduced on -5% over such cities as Achinsk, Novosibirsk, Abakan, Kemerovo (are shown in Fig. 6a) and Krasnoyarsk (outside the map of Fig. 6a), compared to the baseline period. Moreover, the NO 2 reductions have continued over these Siberian cities (reaching − 14%) in 2015-2019, compared to the baseline period, whereas the more recent reduction of NO 2 (− 4%) was registered over Barnaul (Fig. 6b). At the same time, the dynamics over the Kuzbass basin (and related coal quarries), except a minor regional cluster, were strikingly opposite. There was a prominent increase of NO 2 over Kuzbass, whereas the average NO 2 over coal quarries within Kuzbass increased on ~ 7 ± 5% and ~ 15 ± 8% in 2010-2014 and 2015-2019, respectively, compared to the baseline period. The only exception was a minor regional cluster in Kuzbass near Mezhdurechensk, where NO 2 remained nearly unchanged in 2009-2019, compared with the baseline period (0-1% change). We suggest this was due to its topography, as the city is at the highest point in Kuzbass (532 m asl; Fig. S1), which is otherwise mostly flat (< 400 m asl). Another factor might be the lack of coal production plants within Mezhdurechensk, whereas Belovo and Kiselevsk each have two, and Prokopyevsk has one (Table S4 in the Supplementary Material). The latter factor potentially indicates that the combination of coal mining activities (i.e., mines) and coal production facilities (i.e., coal preparation plants) exacerbates the atmospheric NO 2 anomaly in the Kuzbass cities with both coal mining and processing activities.
Novokuznetsk, which is also within Kuzbass, exhibited a slightly different pattern from the increasing NO 2 seen over Kuzbass (Fig. 6). In particular, NO 2 levels fell in 2010-2014 (by 4% relative to 2005-2009), and subsequently returned to the baseline level in 2015-2019. The different pattern for Novokuznetsk is reasonable, as the city is located within Kuzbass but has no coal mines (see Fig. 2) because it specializes in industrial production 28 , where ~ 62% of city enterprises produce, supply or support the production of metal 35 . Interestingly, although Novokuznetsk is not a city with coal-oriented economy, coal is actively used as input in its production of metal as conversion of most or all metal ores to usable metal is highly energy intensive. Metal production facilities use coal to provide energy and the metals are being produced by conversion of coal to coke, where both processes emit substantial NO 2 emissions in the atmosphere. Although the data on metal production of Novokuznetsk is scanty, we analyzed Novokuznetsk inventory-based NO 2 emissions from a previous study 36 . Notably, NO 2 emissions from metal production of West-Siberian Metal Plant (WSMP) account for 84.2% of all the gaseous pollutants of Novokuznetsk 36 in 2014-2018. Most importantly, we found high correlation (r = 0.76) between inventory-based NO 2 emissions from WSMP 36 and our annual NO 2 tropospheric estimates from OMI and, where even higher correlation (r = 0.84) was discerned between WSMP NO2 emissions and coal production from Fig. 4. These www.nature.com/scientificreports/ findings clearly indicate that metal manufacturing of Novokuznetsk is based on coke-intensive input, thereby pointing on a key role of coal in these emissions.

Discussion
The sensitivity of tropospheric NO 2 (measured by OMI) to the mining, production, and transportation of coal in Russia's largest coal basin (Kuzbass) was demonstrated for the first time. A major long-term tropospheric NO 2 anomaly was revealed over Kuzbass (3.22 ± 0.52 × 10 15 molecule/cm 2 ) in the period 2005-2018, indicating substantial gaseous pollution over the region. The anomaly was attributed to the Kuzbass coal basin, based on moderate agreement identified between (1) the spatial distributions of NO 2 and night lights originating from the cluster of coal mines and cites, as well as by (2) the correlation between inter-annual coal production and annual NO 2 levels in Kuzbass (r ≥ 0.60). Unlike the global trend of NO 2 reduction over urban areas (including 25°. This analysis was conducted for retrieving NO 2 over areas with elevated NO 2 concentrations, including cities and the cluster of coal mines in Kuzbass, whereas the hinterlands were less important for this analysis given the lower tropospheric NO 2 over these areas (Fig. 3). To this end, the NO 2 data from areas with relatively low NO 2 were masked out (white areas). Areas with low NO 2 were defined as follows: the NO 2 concentration for a 0. Production in the two latter periods increased by 11% and 30%, respectively, compared with the baseline period. Remarkably, such increases seem to be strikingly proportional to our reported NO 2 increments over the Kuzbass basin during the same period. Assuming a proportional relationship between NO 2 and coal production, a ~ 1.0% increase in coal production is likely to cause a ~ 0.5%-0.6% increase in the NO 2 concentration in the region, where some portion of this increase may have been offset by regulations limiting emissions.
From regional perspective, the demonstrated association between tropospheric NO 2 with coal mining and processing activities over the largest coal basin in Russia is valuable, given the limited opportunity to otherwise assess the environmental impacts of coal-related activities in developing/emerging economies such as Russia, where official information is often inaccurate and atmospheric observations are scarce. These first estimates of substantial tropospheric NO 2 increases over Kuzbass can encourage national and regional policy makers to formulate new pollution mitigation strategies to assess and to minimize the local population's exposure to the adverse cardiovascular and respiratory effects of tropospheric NO 2 and also O 3 , as NO 2 is a precursor of O 3 . Moreover, elevated atmospheric NO 2 may cause indirect adverse environmental effects such as nitrogen enrichment of water bodies via deposition, which compromises the safety of drinking water as well 37 .
Our findings are novel at a global level, as most existing spaceborne studies focusing on coal mining have reported increases in atmospheric CH 4 6,7,38,39 . Although we demonstrated that spaceborne observations of NO 2 can be utilized to attribute NO 2 pollution to previously unreported coal mining and processing over large exploited coal basins, the hypothesis that coal mining itself can be a major source of NO 2 emissions should be evaluated. The tropospheric NO 2 anomaly over Kuzbass could be driven both by direct nitrogenous outgassing from coal mining (e.g., scraping and fracturing) and indirect nitrogenous emissions from coal mining and processing. Such indirect nitrogenous emissions originate from the inefficient compression-ignition diesel engines of coal mining equipment (excavators), coal transportation vehicles (haul trucks) and from nearby facilities for processing coal. Notably, the combustion of heavy fuel by the haul trucks deployed in Kuzbass has already been questioned from an ecological viewpoint 19 . Although new heavy machinery meeting Euro-3 and Euro-4 standards is being delivered to Kuzbass, there remains an abundance of inefficient diesel-fueled equipment prone to emitting NOx 19 . Moreover, as preliminarily indicated, such indirect coal emissions can stem from manufacturing of steel, which uses coke in the production, thereby generating various indirect air pollutants including NO 2 .
We encourage future studies to use Tropospheric Monitoring Instrument (TROPOMI) observations for the high-resolution analysis of NO 2 over Kuzbass as the region is virtually unknown by the TROPOMI research community. Such studies could accurately estimate emissions from coal-related sources in Kuzbass by implementing inversions into atmospheric transport models with high-resolution spaceborne NO 2 observations as the input. In this context, the combination of TROPOMI observations with the emission estimates from inventories can (a) disentangle direct emissions from coal mining and processing from indirect emissions originating from the inefficient fuel combustion of coal mining machinery and coke production of steel industry. Moreover, in this way, one can (b) elucidate the contributions of other human activities unrelated to coal mining to the identified NO 2 anomaly over Kuzbass. Once these aspects are clarified, further research can estimate the coal mininginduced direct nitrogenous outgassing from Kuzbass.

Methods
Proxy data: anthropogenic activities and coal production in Kuzbass. As no institution has consistently reported statistics on coal production in Kuzbass, we collected fragmented data in coal production from various sources including the Administration of Kemerovo Region, the Department of Coal Mining/Production of Kemerovo Region, and the journal Coal of Kuzbass. We compiled coal production statistics for the available years (2006-2018) in Table 1. For further information, see the Supplementary Information and the online compilation of references (https:// gitlab. com/ labzo vskii/ kuzbc oal/).

OMI NO 2 .
OMI is a spectrometer onboard the NASA's EOS-Aura satellite, which is operated in the ultraviolet-visible range (visible channel is used for NO 2 retrieval) since 15 July 2004 at sun-synchronous orbit 40 . We used the retrievals of NO2 tropospheric column, shortly referred as NO2 (Level 3 data, 30% cloud screened product called 'OMNO2dv003') in 2005-2018 period. The number of molecules of NO 2 in atmospheric column were applied with only near-clear sky conditions (cloud radiance fraction < 30%). The spatial resolution of NO2 estimates is 0.25° × 0.25°, produced by the averaged observations (originally ranging from ~ 13 km × 24 km in near nadir to ~ 24 km × 160 km for the observations at the edge of a swath) over a grid cell of interest 41 . The version 3 retrieval stands out with high accuracy of retrievals of NO2 slant column density 42 . This study used the long-term averages (2005-2018) of NO 2 as representative long-term estimates over a broad region, following previous studies 16 and prior indications about the link between coal mining and processing activities and NO 2 emissions 43 . The period of 2005-2018 was selected for OMI based on the availability of coal production data, namely, from the first year with full year OMI NO2 coverage until 2018, the last year for which the coal production in Kuzbass was available. The analysis was extended on one year (2005-2019) for Fig. 6 as we did not need to use coal production as a reference the analysis, shown on this figure. We also calculated annual averages of NO 2 over the region of interest (Kemerovo Region and Kuzbass), and provided the corresponding standard deviations. The data were accessed (on 07.12.2021) using the Giovanni tool of the NASA EarthData service (https:// disc. gsfc. nasa. gov/ infor mation/ tools? title= Giova nni). www.nature.com/scientificreports/ Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.